function [f,df]=f_ola_imp(y,p)
% ejecutar [f,df]=f_ola_imp(x,p)
% donde p=[1,0] calcula f
%       p=[0,1] calcula df
%       p=[1,1] calcula f y df
f=[];
df=[];
if nargin<2
    p=[1,0];
end
global x

if p(1)==1
   f=y-exp(-(1-x*y).^2);
end
if p(2)==1
   df=1-2*x*(1-x*y).*exp(-(1-x*y).^2);
end

